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The two-layer square lattice quantum antiferromagnet with spins i shows a zero-field magnetic 
order-disorder transition at a critical ratio of the inter-plane to intra-plane couplings. Adding a 
uniform magnetic field tunes the system to canted antiferromagnetism and eventually to a fully 
polarized state; similar behavior occurs for ferromagnetic intra-plane coupling. Based on a bond 
operator spin representation, we propose an approximate ground state wavefunction which consis- 
tently covers all phases by means of a unitary transformation. The excitations can be efficiently 
described as independent bosons; in the antiferromagnetic phase these reduce to the well-known 
spin waves, whereas they describe gapped spin-1 excitations in the singlet phase. We compute the 
spectra of these excitations as well as the magnetizations throughout the whole phase diagram. 



I. INTRODUCTION 

Bilayer quantum magnets have attracted much inter- 
est in experiment and theory in recent years, especially in 
the context of quantum Hall systems and of quasi- two di- 
mensional transition metal oxides. On the experimental 
side, quantum Hall systems are especially suitable for 
investigating zero temperature quantum transitions be- 
tween states with different spin magnetizations. In par- 
ticular, bilayer quantum Hall systems at filling fraction 
V = 2 d"^] can be tuned between a fully polarized, 
ferromagnetic state and a spin singlet ground state as 
function of the layer distance. It is now well established 
ll|iil"0 that there is not a first-order ferromagnet-singlet 
transition, but rather an intermediate phase with canted 
spin ordering bounded by second-order transitions. The 
spin degrees of freedom therefore appear to be well de- 
scribed by a bilayer quantum spin model On the 
other hand, transition metal oxides like cuprates are 
known to form two-dimensional structures where the low- 
energy spin dynamics is well described by a Heisenberg 
model. These materials consist either of a single plane 
or a stack of copper oxide planes with intervening charge 
reservoir layers. Prominent examples for bilayer com- 
pounds are YBa2Cu307_a; and Bi2Sr2CaCu208+5 which 
show a number of unusual properties Note that in 
the high-temperature superconductors the interlayer cou- 
pling is relatively weak compared to intralayer processes; 
undoped bilayer systems like YBa2Cu306 are not driven 
into a singlet ground state. However, the latter is also 
possible: recently the material BaCuSi2 06 containing 
strongly coupled bilayers has been discovered and inves- 
tigated [0 , it shows a spin gap due to strong antiferro- 
magnetic exchange interaction between the two layers. 

On the theoretical side, the bilayer Heisenberg mag- 
net has attracted a lot of interest [^,[ll|-|25t because it is 
a simple model for studying the interplay of long-range 



magnetic order and quantum disorder. Quantum-critical 
behavior associated with such a magnetic instability has 
been observed, e.g., in the cuprate superconductors over 
a wide range of doping levels and temperatures ||2^ . 

We start by describing the bilayer quantum Heisenberg 
model. The system consists of two planes of nearest- 
neighbor S — ^ Heisenberg models with coupling con- 
stant J|| which can be either antiferromagnetic ( Jy > 0) 
or ferromagnetic (J|| < 0). The spins of corresponding 
sites of each layer are coupled antiferromagnetically with 
a coupling constant Jj^ > 0. The Hamiltonian reads 

^ = ^ Sii -8^2 + J^ll ^ Sim-Sjm — Sim (1) 

i {ij)rn im 

where Sim are the electronic spin operators. The in- 
dex i denotes rungs of the bilayer lattice, and m — 1,2 
labels the planes. The summation (ij) runs over pairs 
of nearest-neighbor rungs. The external magnetic field 
B is homogeneous, the Zeeman interaction factors gj^s 
(gyromagnetic ratio and Bohr magneton) have been ab- 
sorbed in the definition of B. 

In the following, we briefly discuss the main features of 
the ground-state phase diagram of the model (^, which 
has been given in two recent papers ||8|,|l5|| and is shown 
in Fig. |l]. In the absence of a magnetic field the bilayer 
Heisenberg antiferromagnet is known to exhibit quantum 
phase transitions between a disordered singlet phase and 
long-range ordered (LRO) phases. In the limit of strong 
interplanar coupling, J± ^ |J|||, the system consists of 
weakly interacting rung singlets, and the ground state 
possesses the full symmetry of the Hamiltonian. The 
spin excitations are triplet modes with a minimum en- 
ergy gap A, there is no magnetic LRO. In the oppo- 
site case of large |J|||/Jj_, the system possesses LRO at 
T = 0, the ordering wavevector is (0, 0, tt) [(tt, tt, tt)] for 
ferromagnetic [antiferromagnetic] J||. The SU(2) symme- 
try of H is broken, the low-lying excitations are doubly- 
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FIG. 1. Ground state phase diagram of 7i (jj) as deter- 
mined from series expansion studies, see Refs 15po| . The 
arrows denote the spin orientations, the in-plane ordering 
wavevectors are given in brackets. The dashed lines are the 
boundaries of the singlet phase obtained from the lowest-order 
boson approximation as described in Sec [II, 



degenerate Goldstone spin waves. The quantum tran- 
sitions between the singlet and the two ordered phases 
are of the 0(3) universaHty class ||ri|,|l^,|l8|,|3| and occur 
at critical ratios J||)ci,2- For the antifcrromagnetic 
case, quantum Monte Carlo calculations [pT|jl^j25|] , se- 
ries expansions |p^ , and the diagrammatic Brueckncr 
approach yield an order-disorder transition point 
of (J||/Jl)c1 = 0.396. Bond-operator mean-field the- 
ory applied to the bilayer Heisenberg AF gives 
a transition point of (J||/J_l)c1 = 0.435. Note that 
Schwinger boson mean-field theory [T^ ] predicts a value 
of (J||/Jx)ci ~ 0.22, and also self-consistent spin-wave 
theory [||,|9|, which yields (J||/Jj.)ci « 0.23, fails to re- 
produce the numerical results. As Chubukov and Morr 
I p^ have pointed out this discrepancy is due to the ne- 
glect of longitudinal spin fluctuations in the conventional 
spin-wave approach. For ferromagnetic J|| < 0, the series 
expansion result is ( Jy / J±)c2 = —0.435 which agrees well 
with the value from bond-operator mean-field theory. 

For large enough external field, B = Be^, the exact 
ground state is the fully polarized ferromagnetic (FPF) 
state, i.e., the state with all spins pointing in the direc- 
tion of the applied field. The first excited state is a single 
spin flip and its excitation energy can be determined ex- 
actly: Wk = B — J± — J\\{2—coskx—cosky). For ferromag- 
netic [antifcrromagnetic] intra-plane coupling J[| the min- 
imum of Wk is located at in-plane wavevector Q = (0, 0) 
[Q = (tt, tt)]. The stability boundary of the FPF state is 
given by the point where the minimum excitation energy 
vanishes, this yields the exact expression for the upper 
critical field, 3^2 = [Bc2 — 4J|| + J±] for ferromag- 
netic [antifcrromagnetic] J|| . Single spin flips condense 
at this boundary, leading to canted spin ordering. The 



transition is of second order and in the universality class 
of the dilute Bose gas quantum-critical point with z = 2 
[]6|-p|, p"5p^ ] . For intermediate magnetic flelds a canted 
spin phase is established which breaks the rotational sym- 
metry of the Hamiltonian about the z axis, leading to one 
linear dispersing Goldstone mode corresponding to a ro- 
tation of the order parameter in the x-y plane. Both the 
expectation value of the uniform magnetization in field 
direction and the Q-staggered magnetization perpendic- 
ular to the field are non-zero. 

For small uniform fields, B < A, the disordered spin 
singlet ground state discussed above remains unchanged 
due to its finite excitation gap A. The effect of the mag- 
netic field is simply to split the degenerate triplet exci- 
tations due to Zeeman coupling, Wk ^ t^k ^ mB, where 
m — —1,0, 1 is the Sz quantum number. The stability 
boundary of the singlet phase is determined by the van- 
ishing of the excitation gap; the corresponding transition 
is in the same universality class as the boundary of the 
fully polarized phase, also leading to the same canted 
phase with ordering vector Q. 

The purpose of the present paper is to propose a trans- 
formation of the usual bond operator basis which allows 
an efficient description of all phases described above. In 
the simplest approximation, the excitations can be de- 
scribed as independent bosons. In the limit of weak inter- 
plane coupling, the approach then reduces to linear spin- 
wave theory, with the surplus of correctly describing the 
gapped longitudinal amplitude mode which corresponds 
to inter-plane quantum fluctuations and becomes impor- 
tant near the transition to a singlet state. The paper is 
organized as follows: In Sec || we introduce the bond op- 
erator transformation which is employed to describe the 
magnetic phases of the bilayer model. Sec [II is devoted 
to the boson approximation which can be used to diago- 
nalize the Hamiltonian via a Bogoliubov transformation. 
Results for ground state magnetizations and excitation 

A summary and discussion 



spectra are given in Sec IV 



of possible improvements beyond the independent boson 
approximation close the paper. 



II. GENERALIZED BOND OPERATORS 

The four spin states per rung can be conveniently de- 
scribed via a bond operator representation j2^ of the two 
spins Sii and 8^2 of each rung i. We introduce bosonic 
bond operators for creation of a singlet and three triplet 
states out of the vacuum |0): 

4\0) =^(4,T4,i-4,i4T)|0) = ^(lti>-UT)) (2) 

4io) =^(4_,4T-4,,4i)io) = ^(itT)-Ui)) 
4 io) =-!=(4^^4^^+4_^4j |o) = -j= (iTT)+Ui)) 
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where c|„j^ are creation operators for electrons at site 
im with spin a. In terms of the bond operators the spin 
operators can be expressed as 

A local constraint of the form 

4s^+Y.^lUc. = l (4) 

a 

has to be imposed on each rung to ensure that the physi- 
cal states are either singlets or triplets. The ground state 
for vanishing intra-plane coupling, J|| — 0, is given by the 
product state of singlet bonds on each rung, i.e., 



(5) 



In this limit, the excitations are localized triplets with an 
energy gap J±. 

Substituting the bond operator representation of the 
spin operators into the original Heisenberg model we 
obtain the following Hamiltonian pT 



i ia 

^1 = Y E i^c^^i-s, + h.c.) 



(6) 



In the following the basic idea of describing the vari- 
ous magnetic phases will be discussed. Starting from the 
singlet phase, the magnetically ordered phases can be de- 
scribed by different triplet boson condensates . 
The FPF phase with polarization in z direction requires 
a condensate of {t^ + ity) bosons; a Neel state with 
staggered magnetization perpendicular to the field cor- 
responds to a condensate of or ty on top of the sin- 
glet state. Neglecting inter-rung fluctuations, an ansatz 
wavefunction can be written as 

l^o) ^ exp(z//^4i„)exp(A^e'«^"4sO|</'o) 

i i 

= n(4 + Ae^'^^-(4+^Kl))|0). (7) 



The real parameters A and ^ are the condensation ampli- 
tudes, Q denotes the in-plane ordering wavevector. This 
ansatz contains two subsequent unitary transformations 
of the singlet product state |(/)o) the A term creates a 
spin-density- wave condensate with ordering vector Q and 
quantization axis in x direction [thus explicitely break- 
ing the U(l) symmetry of the z axis rotation], the /i term 
introduces canting and a non-zero z magnetization into 
this state. For = and Q = (vr, tt) which corresponds 
to Jj| > and i? = 0, this wavefunction reduces to a state 
interpolating between the singlet and the Neel state; it 
has been recently used to describe the dynamics of holes 
doped into a bilayer system A finite external field 
B larger than the spin gap will lead to a finite uniform 
magnetization described by finite /i; for large B the sys- 
tem is driven into a fully polarized state with A — > cxd 
and ^ = 1. The intervening canted phase will be charac- 
terized by finite, non-zero A and < /i < 1. 

A related product state has recently been employed for 
the description of bilayer quantum Hall systems ||^ , with 
the inclusion of charge fluctuations and additional dis- 
order, but without inter-rung spin fluctuations. We will 
show here that such spatial correlations can be efficiently 
described in a simple harmonic theory. 

For a proper description of fluctuations around the 
product state |0o) (0) it is convenient to transform the 
basis states on each rung. We replace the basis operators 
{si,tia} by 
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(8) 



The constraint takes the form sjsi + ^la^ia = 1- 
The new product state \(f)o) can now be written as 



i0o)-n-?iio) 



(9) 



which reduces to the singlet product state for A = 0. Ex- 
citations out of this state are given by il^Si. At this point 
it is worth emphasizing that the use of a single conden- 
sate (of "mixed" s bosons) and the basis transformation 
to excitation operators being orthogonal to s is crucial 
for a correct description of the ordered phases. 

We can now insert the new basis operators into and 
obtain the expression for the transformed Hamiltonian 
7Y: 
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ia i 
i 

i 

+ 5^ iJ(ss; ss) sJ4^*^J' 

+ 5] ss) t^ia^lj^Sj + h.c. 

{ij)a 

+ Y H{xy; ss) tl^t'^ySiS^ + h.c. + j 



-ff(as; sa) il^ijas''jSi + h.c. 

{ij)a 

iJ(a;s; sy) tl,^ijys'jS.i + h.c. + i ^ j 

i?(a;s; xs) iL^iajS^Sj + i ^ j 
51 ^(y^'^ y^) tlyi^ySpj + j 

(ij) 
(ij) 



+ Y, ^i^s'^ ss) (i|a;Sisjsj + s]ii2:s]sj^ + i ^ j 
+ Y ^^y^'' + ^l^w^pj) + j 

+ Hrcst . (10) 

The coefficients •) depend on tlie basis parameters 
A and ^ as well as on J||, J_l, and i3; explicit expression 
are given in Appendix A. Titcst contains more compli- 
cated terms with three and four excitations. In the case 
of = the Hamiltonian reduces to the Hamiltonian 
derived for the zero- field case [|T|. It can be seen that 
Til contains creation, hopping, and conversion terms of 
the three types of excitations {tix,tiy,tiz}. Contrary to 
the original Hamiltonian Tii where the triplet excitations 
can only occur in pairs, here also single tix and tiy ex- 
citations can be created. The effect of these terms is 
directly related to the basis parameters A and /i and will 
be used to determine their values. In contrast to the case 
of zero field where only pairs of the same type of excita- 
tion are created out of |0o), for any finite B also pairs of 



the structure t^iJ^jySiSj are present. 

The exact but lengthy representation (^0|) can now 
be used for a variety of approximation schemes. The 
most important point is to choose values for A and /i 
in such a way that the product state |0o) @j is a rea- 
sonable starting point for approximations. Fluctuations 
are described by the hard-core bosons {tix,tiy,tiz}, and 
a "good" choice for \(j)o) should ensure that the ground- 
state density of the bosons is small. The lowest-order 
approximation corresponds to independent bosons with 
a bilinear Hamiltonian, it is formally given by ignoring 
both the hard-core constraint as well as Tirest , and as- 
suming a complete condensation of the generalized "sin- 
glet", i.e., (s) = s = 1. This approximation will be 
discussed in the following sections, it is shown to corre- 
spond to usual linear spin-wave theory in the decoupled 
plane limit with B = 0. Possible improvements include 
(i) taking into account the site-averaged constraint by 
introduction of a chemical potential Aq, and treating Aq 
and s as variational parameters in the spirit of bond- 
operator mean-field theory (ii) a mcan-field-like 
factorization of the higher-order boson interaction terms 
(iii) a diagrammatic treatment of the hard-core boson in- 
teraction in the framework of Brueckner theory [p^ . In- 
dependently, methods focussing on higher-order local ex- 
citations based on cumulants |31| or the coupled-cluster 
technique [ p3[ known from quantum chemistry may be 
applied to po|). 



III. INDEPENDENT BOSON APPROXIMATION 

In this and the following section, we will discuss 
the lowest-order boson approximation for the Hamilto- 
nian ( p^) . For this purpose we treat the excitations 
{Ux-iUyiUz} as independent bosons, i.e., neglect the 
constraint which restricts the Hilbert space per rung. 
The (generalized) singlet is assumed to fully condense. 



(s) = 1. It is known from the zero-field case [|_4 31 
that this approximation is controlled by the existence of 
a small parameter, namely the density of (generalized) 
triplet excitations {t[Jtia). With this in mind, all terms 
of the Hamiltonian containing more than two excitation 
operators tia will be neglected. 

The parameters of the unitary transformation A and /i 
have not yet been determined. Within the approximation 
of independent bosons these parameters are chosen so 
that the prefactors of the terms creating single Ua bosons 
out of |(/)o) vanish. The physical meaning of this is easily 
understood: terms creating single bosons would change 
the condensate densities (and therefore alter the effective 
values of A and /i), however, our aim is to fully account 
for the boson condensation by the transformed basis state 
1 00) • One arrives at the following non-linear equations for 
the basis parameters A and 
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[±4J|| - J±{l+n^) + 2Bfi\{l+fi^) 
[ ± 4J,| + Jj_(l + Ai2) - 2Bfi] (1 + ^2) + 8J||Ai2 



(11) 



= B(1-/)(1 + A2) =f 4J||^(1 + ^2) _ 4J||AV(1-Ai^) 

where upper (lower) sign refers to J|| > (J|| < 0). No- 
tably, these equations also follow from a variational prin- 
ciple, i.e., they are obtained by minimizing (0o|^|</'o) 
with respect to A and fi. This clarifies the meaning of 
|(/)o) as best variational state without inter-rung fluctua- 
tions. 

The system of nonlinear equations (|l^) allows for the 
computation of the phase boundaries. The rotational 
invariant spin singlet phase is characterized by the value 
A = 0. Solutions for fi are then obtained for B < Bd, 
where Bd is the stability boundary of the singlet phase: 



Bel = - 4J|| ( J|| > 0) , 
J± + 4 J|| ( J|| < 0) . 



Bel — 



(12) 



It can be seen that the singlet phase exists for | J||/J^| < 
0.25. At the critical ratio {J\\/J±)c = ±0.25 the char- 
acter of the zero-field ground state changes from spin 
singlet to long-range order. As discussed in Ref |23|, this 
quasi- classical value of {J\\/J±)c is necessarily smaller 
than the exact one (quantum fluctuations stabilize the 
inter-plane singlet phase); as noted in the introduction 
the critical values obtained by large-scale numerics are 
{J\\/J±)ci = 0.396 and (J||/Jj.)c2 = -0.435. The in- 
clusion of higher order terms beyond the present linear 
"spin-wave" approximation leads to a better agreement 
with the numerical values for the critical coupling ratio 

ill- 

For values of the field with Bd < B < Bc2 the ground 
state is a canted phase with < /i < 1. The phase 
boundary to the ferromagnetic phase is obtained when 
the limit A — > cxd and /i ^ 1 is taken. We arrive at an 
upper critical field of 

Bc2 = J±+ 4 J|| ( J|| > 0) , 

Bc2 = J± {J\\ < 0) (13) 

which are the exact results. 

We can now proceed with analyzing the fluctuations 
around \(t>o)- With the approximations described above, 
we arrive at a bilinear Hamiltonian: 



C'k [tl jky - ilyii^x 



momentum k. The expressions for the coefficients 
^k, • • ■ , Dk can be found in Appendix B. The Hamilto- 
nian (|l^ can be easily diagonalized with a Bogoliubov 
transformation, leading to new generalized excitations 
denoted by Tkc, a — x,y,z. The ground state j'i/'o) is 
simply the vacuum of the bosons Tka- The excitation 
energies of these bosons are given by 

V 2 



-^kx^ty ^ (Ck ^ ^k) 

+2 ((^kx^kiy — B]^xBky) (Ck + i?k) 

— 2 (Akrc-Sky — AkyBi^x) CkDi^) 



^kz — ^kz B^^ , 



(15) 



the terms 



■— ^kx 



BZ^ and uoT, 



^ky 



Biy are 



— — - 

''kx ■ '^kx ^kx '^^^^ ky 

simply the eigenenergies of the zero field limit. In zero 
field, we find either three degenerate gapped modes in 
the singlet phase, or two degenerate transverse acous- 
tic spin-wave modes and a gapped spin-amplitude mode 
in the long-range ordered phases. A finite external field 
lifts these degeneracies; the canted phase supports a sin- 
gle Goldstone mode due to the broken U(l) symmetry. 
Explicit results for the dispersion relations are shown in 
the next section. 

Finally the ground state energy per site is given by the 
expression 



Ey = 



2N 



^^" + ^E("ka-Aka) (16) 
ka 



(14) 



where Eq is the energy per site of the product state |(/)o), 
Eq = (^o|H|0o)/2iV. 

The described harmonic approximation preserves 
SU(2) symmetry in the singlet phase; in the ordered 
phases it can be considered as linear spin-wave theory 
for the bilayer problem with inclusion of longitudinal spin 
fluctuations p8| , ^j34| - these are crucial for describing 
the magnetic properties near the boundary to the singlet 
phase. For zero field and vanishing inter-plane coupling, 
B = — Q, the results of the present approach are 
equivalent to the ones of linear spin-wave theory. For 
instance, in the antiferromagnetic case, the magnetiza- 
tion takes 60 % of its classical value. There, the gapped 
amplitude mode is dispersionless and does not influence 
the ground state properties, and its spectral weight in 
neutron scattering vanishes. 



Here, ika are the modified basis operators which have 
been Fourier transformed with respect to the in-plane 
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Magnetic field B [J^+ 47||] 

FIG. 2. Uniform (top) and staggered magnetization (bot- 
tom) as function of the external magnetic field B/Bc2- The 
ratio of the coupling constants J\\/J± ~ 0.2 which places 
the zero-field system into the interlayer singlet phase. Solid: 
present approximation (independent bosons). Dashed: prod- 
uct wavef unction |0o) only. 

IV. GROUND STATE PROPERTIES AND SPIN 
EXCITATIONS 

In this section we present the results obtained from 
the independent boson approximation. We discuss the 
staggered and uniform magnetizations and the dispersion 
relations as functions of the external magnetic field in 
order to characterize the various magnetic phases. For 
this purpose we consider the system for small intra-plane 
coupling |Jj[|, i.e., in the singlet phase, as well as for 
weakly coupled planes (small J±). 

A. Magnetization 

We start with the uniform and staggered magnetiza- 
tions which are obtained as expection values of the cor- 
responding spin operators 

i 
i 

where Q = (0, 0) [Q = (tt, tt)] for ferromagnetic [antifer- 
romagnetic] intra-plane coupling J|| as above. The stag- 




0.2 0.4 0.6 0.8 1 
Magnetic field B [/J 

FIG. 3. As Fig. ^ but in the singlet phase with strongly 
coupled ferromagnetic planes with J\\/J± ~ —0.2. 

gered magnetization in the zero-field case has been calcu- 
lated earlier (see Fig 2 of Ref ||3|]); it shows a good over- 
all agreement with the series expansion data of Ref jl^] 
if plotted as function of {J\\ / J±) / {J\\ / J±)c - the critical 
coupling ratio itself is underestimated in the independent 
boson approach as discussed above. 

Turning to the finite-i? case, we show uniform and 
staggered magnetization in Figs Both are zero in 

the singlet phase, B < Bd. Varying the magnetic field 
from the lower critical value B^i to the upper critical 
field Bc2 the uniform magnetization increases from and 
saturates at Bc2 , whereas the staggered magnetization is 
non-zero only in the canted phase, Bci < B < Bc2, pass- 
ing a maximum in between (Figs H, Near the critical 
fields for B > B^ or B < Bc2, we have M„ ^ \B ~ Bci,2| 
andM,t~ |5-B,i,2r/2. 

For small J^, i.e., in the zero-field ordered regime, the 
staggered magnetization M^t is finite B — Q but is 
suppressed if B is increased up to the critical field i3c2, 
see Fig^. Interestingly, Mst increases (!) at small fields - 
here the suppression of quantum fluctuations due to the 
field [e.g. one Goldstone mode becomes gapped) has a 
stronger effect than the change in the canting angle. The 
uniform magnetization again increases continuously with 
B until it saturates at Bc2- For illustration, we also show 
data for J\\/Ji. — —0.25 which places the model at the 
zero-field critical point between the singlet and Q = (0, 0) 
Neel phases. Here, M„ increases with B^ whereas M^t 
raises linearly with the applied field, in agreement with 
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0.2 0.4 0.6 0.8 1 
Magnetic field B [J^+ 4J^^] 

FIG. 4. As Fig. ^, but for weakly coupled antiferromag- 
netic planes with J\\/J± = 10. 




Magnetic field B [J J 

FIG. 5. As Fig. ^, now at the zero-field critical point be- 
tween the singlet and "ferromagnetic" Q — (0, 0) Neel phase, 
J||/Jx = -0.25 (see Fig. §). 



the results of Ref Q . 

In Figs ||-^, we have also shown the magnetization val- 
ues calculated with the product state |<^o) only (dashed 
lines) - this is equivalent to neglecting intra-plane (inter- 
rung) fluctuations altogether. It can be seen that the 
effect of these spatial quantum fluctuations is largest for 
antifcrromagnetic Jy and large staggered magnetization, 
this includes the case of the single-layer antiferromagnet 
where quantum fluctuations are known to reduce Mgt 
by 40%. In contrast, in both the ferromagnetic regime 
(where quantum fluctuations are weak in general) and 
near the spin singlet phase (where inter-plane quantum 
fluctuations dominate), the effect of intra-plane fluctua- 
tions is weaker. 

All phase transitions found are of second order, as ex- 
pected, except for the special point Jy = 0, B = J±, 
where the singlet-FPF transition is a simple level cross- 
ing for independent rungs. Due to the mean-field char- 
acter of the present approximation, all critical expo- 
nents have mean-field values. (It is easily seen that 
the present harmonic approximation becomes exact in 
the limit of large in-plane coordination number, see also 
Ref [^.) For the actual two-dimensional model, these 
mean-field exponents are incorrect for the zero-field tran- 
sitions which have z = 1 and therefore obey the non- 
trivial exponents of a three-dimensional classical Heisen- 
berg model, but they apply (ignoring logarithmic correc- 
tions) to the field-induced transitions from and to the 
canted phase which are at the upper-critical dimension 
{z = 2) 

We note that the phase boundaries obtained in Sec III 
using the product state \(t>o) are not changed by the inclu- 
sion of inter-rung fiuctuations at the lowest boson level 
employed here. It is, however, clear that boson interac- 
tions change this picture fl^ j. 



B. Magnetic excitations 

We continue with the discussion of the dispersion re- 
lations of the magnetic excitations, plotted in Figs 
for a number of values of the external field. In Fig ^ we 
show data for strongly coupled planes and antifcrromag- 
netic Jy . In the zero-field limit (Fig ^) the excitations 
are threefold degenerate, gapped, and have a dispersion 
minimum at Q = (tt, tt). At finite B, the degeneracy is 
lifted due to the Zeeman coupling. At the lower criti- 
cal field, B = Bel (Fig |^b), the excitation gap closes, 
and there is a single critical mode with quadratic disper- 
sion. In the canted phase, Bd < B < Bc2 (not shown), 
there is a single Goldstone mode with a linear dispersion 
around Q = {tt,tt), corresponding to the broken U{1) 
symmetry. At the phase boundary to the ferromagnetic 
(FPF) phase, B = Bc2 (Fig ||c), the lowest mode be- 
comes again critical with a quadratic dispersion (note 
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(0,0) (71,0) (71,31) 

In-plane wave vector k 



(0,0) 



FIG. 6. Dispersion relations of the magnetic excitations 
for strongly coupled antiferromagnetic planes, J\\/J_l = 0.2. 
The three panels show the three modes for different magnetic 
fields B = BjBci = B/(4J|| + J±), (a) B = 0, (b) B = 0.44 
(here B = B^i), [c) B = I [B = B,2). 
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FIG. 8. As in Fig. ^, here for weakly coupled antiferro- 
magnetic planes, J\\/ Ji_ ~ 10. 
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(71,0) (7l,7l) 
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FIG. 7. As in Fig. ^ but for ferromagnetic in-plane cou- 
pling, J\\/J± ~ —0.2. Note that the magnetic field values in 
units of Bc2 are again B = 0, 0.44, and 1, but here Bc2 ~ J±, 
i.e., B = B/J±. 
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FIG. 9. As in Fig. ^, finally for the coupling constant ratio 
J\\/J± = —0.25, corresponding to the zero-field critical point 
on the ferromagnetic side of the phase diagram (Fig. Q). 
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z = 2 at both field- induced transitions). As the mag- 
netic field is further increased, a gap opens and the FPF 
phase is stabilized. In the FPF phase, one of the exci- 
tations is dispersionless, corresponding to a longitudinal 
local fluctuation which turns both spins on a rung down. 
In the case of ferromagnetic in-plane coupHng, J|| < 0, 
the picture is similar, but the dispersion minimum is now 
found at Q = (0, 0), see Fig 0. 

The results for weakly coupled antiferromagnetic 
planes are plotted in Fig ||. We start with a discus- 
sion of the zero field case, shown in Fig |^a. Contrary 
to the case of strong inter-plane coupling there is now 
a nearly dispersionless excitation - this is the before- 
mentioned gapped (longitudinal) spin-amplitude mode. 
It corresponds to flipping both spins on a single rung 
and has a quantum number S*^ = 0, consequently, it re- 
mains unchanged (to linear order) when the magnetic 
field is turned on. The two other modes are degenerate 
Goldstone bosons in the zero-field limit with a linear dis- 
persion around the ordering wavevector Q = (tt, tt). For 
decoupled planes, J±_ 0, momenta (0, 0) and (tt, tt) are 
degenerate, and for small J± the gap at (0, 0) is propor- 
tional to J_L. With increasing Jx, the spin-amplitude 
mode acquires a dispersion, until it merges with the 
acoustic spin- wave modes at the zero-field critical point, 
leading to a degenerate triplet of critical modes with lin- 
ear dispersion (not shown). Turning to the case of a finite 
uniform field B, we note that arbitrarily small B estab- 
lishes canted ordering here. For Q < B < Bc2 (Fig ||b), 
the mode degeneracy is again lifted and leaves a sin- 
gle gapless Goldstone excitation with linear dispersion. 
At B = Bc2 (Fig ^), this mode becomes critical with 
quadratic dispersion as discussed above. 

Finally, Fig ^ shows dispersions for the coupling ra- 
tio corresponding to the zero-field ferromagnetic critical 
point, J||/J± = —0.25. For B = (Fig ||a), we have a 
triplet of degenerate critical modes with minimum (con- 
densation) wavevector Q = (0, 0) and a linear dispersion 
reflecting z — 1. With application of a finite field, we 
enter directly the canted phase with a single Goldstone 
mode, so for Figs Hb,c the above discussion applies. 

At this point a remark about the physics of longitudi- 
nal magnon mode is in order. As mentioned, this mode 
is not included in conventional spin-wave theory for the 
antiferromagnet. The reason is that spin- wave theory 
has the character of a l/S* expansion, and longitudinal 
fluctuations are suppressed in the S* — > cx) limit. This 
implies that any theory based on a (finite-order) large- 
S expansion can never treat longitudinal and transverse 
spin fiuctuations on equal footing [ p^ , which is, however, 
crucial near the critical point in the bilayer model where 
longitudinal and transverse fluctuations become indistin- 
guishable. In contrast, our approach preserves the quan- 
tum nature of the 5* = 1/2 spins, but has the spirit of a 
large-0 theory, where z is the in-plane coordination num- 
ber. 



We close this section with some general remarks. The 
lowest-order boson approximation neglects interactions 
between the modes, therefore all excitations including the 
amplitude mode are undamped, and no decay into the 
two-particle continuum is possible. The description of 
such damping processes requires the inclusion of fourth- 
order boson terms, and is deferred to future investigation. 

The low-energy properties of the spin excitations ob- 
tained in the present work are in agreement with hydro- 
dynamical considerations. The spin-wave velocity c in an 
ordered phase, deflned as Oq = c |q — Q| near the order- 
ing wavevector Q, has to fulflU the relation — Ps/x± 
where Ps is the spin stiffness and X-L denotes the uni- 
form susceptibility for a field perpendicular to the or- 
dering direction (here x). At the zero-field transitions, 
J\\/J± = ('^||/'^±)ci,2j both Ps and x± vanish whereas c 
remains non-critical as expected for an 0(3) transition. 
This behavior is correctly reproduced in our calculations 
(Figs H and ^, note x± = dMu/dB. (In contrast, at 
the field-induced transition points the spin-wave velocity 
c vanishes, leading to a quadratically dispersing critical 
mode at i? = i?ci,2-) 

The special properties of the multicritical points at 
B — 0, J\\/J± = {J\\/ J±)ci,2 have been discussed in 
Ref 1^. One consequence of universal scaling arguments 
is that the energy gap of the B = amplitude mode in 
the ordered phase close to the transition is given by the 
spin stiffness ps] this result is special to d = 2 where 
Ps has the dimension of energy. We also note that a 
well-defined amplitude mode requires that the low-energy 
theory for the zero-fleld ordering transition is below its 
upper-critical dimension (here dc,up = 3) so that the in- 
teraction term of the effective 0** theory is relevant in 
the renormalization group sense. Our mean-field theory 
gives an amplitude mode independent of the spatial di- 
mension d, but from the above it is clear that even the 
low-energy part of this mode will be overdamped from 
boson interactions in higher dimensions, d > 3. 

The results obtained with the present method for van- 
ishing inter-plane coupling, J± = 0, are consistent with 
recent calculations using linear spin-wave theory ||35|] for 
a single-plane antiferromagnet. 

V. CONCLUSION 

In this paper, we have proposed a unitary transforma- 
tion of bond-boson operators appropriate for the descrip- 
tion of various magnetic phases, especially with canted 
spin order. The resulting boson Hamiltonian has been 
obtained exactly; it allows for a number of approxima- 
tions. We have applied this method to the bilayer Heisen- 
berg magnet in a uniform external field. Using the lowest- 
order approximation of free bosons we obtained analytic 
expressions for the phase boundaries, the quasiparticle 
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excitations and both uniform and staggered magnetiza- 
tion. The harmonic approximation is controlled by the 
smallness of the density of generalized bosonic triplet ex- 
citations, this parameter has been found to be smaller 
than 0.15 throughout the whole phase diagram. In the 
limit of vanishing inter-plane coupling, the boson ap- 
proximation reduces to linear spin-wave theory; near the 
transition to the singlet phase it has the advantage of 
incorporating longitudinal spin fluctuations. The phase 
diagram (Fig |l|) is adequately described; the locations 
of the zero-field transitions are {J\\/ J±)ci,2 = ±0.25, 
these values are smaller in magnitude than the ones ob- 
tained by more accurate numerical and analytical meth- 
ods [{J\\/J±)ci = 0.396 and (J||/Jj_)c2 = -0.435]. These 
deviations arise from the neglect of boson interactions 
which stabilize the singlet phase. The approach pre- 
sented here preserves all general properties expected from 
hydrodynamics, i.e., it gives the correct number and low- 
energy dispersion of Goldstone modes in the symmetry- 
broken phases. 

The present approximation can be systematically im- 
proved as discussed in Sec Two promising possibili- 
ties are either the inclusion of fourth-order boson terms 
as in non-linear spin-wave theory or the treatment of the 
hard-core boson constraint by introduction of an (infi- 
nite) on-site repulsion which is treated diagrammatically 
using Brueckner theory |Q . (The latter approach yields 
very accurate results for the zero-field bilayer antiferro- 
magnet.) We expect that the inclusion of interactions 
will induce a damping of the high-energy spin excitation 
modes, but will not change the low-energy properties pre- 



sented in Sec [V which are protected by symmetry con- 
siderations and hydrodynamics. Once interactions are 
included, a detailed study of the dynamic structure fac- 
tor as measured in inelastic neutron scattering, including 
the spectral weights of the excitation modes described 
here, would be interesting. 

We close by mentioning possible applications of the 
present approach. As in Ref it can be used for 
studying a hole- or electron-doped system; in this case 
the ground state wavefunction derived here can be used 
as a background state for the carrier motion. Other in- 
teresting prospectives include the study of disorder |^] 
and the application to frustrated systems in the presence 
of an external field. A nice example here is the recently 
discovered triangular bilayer S = 1 system Ba3Mn2 08 
which shows a spin gap, magnetization plateaus, and in- 
teresting frustration effects [Bq| . 



APPENDIX A: HAMILTONIAN IN THE 
MODIFIED BOND BOSON BASIS 

In the following, we list the first terms appearing in 
the transformed bilayer Hamiltonian H ( |Tc| ) for antifer- 
romagnetic J||. 



1 + At V 1 + A 



Hiy;y) = J^ + 2B 



H{x; s) = {Jj_-2B 



l + /i' ' 



H{z-z) = Ji 



1 + 1 + A2 



H{y;s) = B 
H(x-y) 



1 + VTT^ ' 
-B 1-^2 



H{xx] ss) 
H{yy;ss) 
H{zz; ss) 



{1 + X^f {1 + fj.^f 



(1-M^) [(1 + M^)+A^(1-M^)] 
2 (1 + A2) (1 + ^2)2 

(l + /i2)+A2(l-^2^ 
2 (l + A2)(l + ^2) • 



These and the remaining terms can be obtained with the 
help of Mathematica by inserting the transformation (|^) 
into n (|). 



APPENDIX B: BILINEAR HAMILTONIAN 

We list here the coefficients of the bilinear boson 
Hamiltonian (^4|) which enter the expressions for the dis- 
persions given in (p^. For shortness, we restrict our- 
selves to the antiferromagnetic case, Jy > 0, and use the 
abbreviation 7k = 5 {coskx + cosfcj,). 
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1 + 1 + A 



(1 + A2)2 (1+^2)2 



2>/||7k' 



^1-A2)^(1 + 2m^) + (1 + A2)>4 

(l + A2)2(l+;,2)2 
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A-liy — J± 



2Bfi- 



1 + 2A2 



8J| 



1 + A2 ' '^(l + A2)(l+^2) 

A2 (1 - 2X^fi^) 

'(1 + A2)2 (1+^2)2 



2>/||7k 



8Ji 



(1+/.2)^_A2 (1-^2)2 
(1 + A2) (1+^2)2 ' 

1 A2 

1+A2 +^^^(1 + A2)(1+/Z2) 
A2(l+/i2(l_A2)) 



(1 + A2f (l+/.2)2 ' 2^li^ni + A2; ' 



Ska; = 2J[|7k- 



fl _ XA'^ _ n ^ \2^2 ,,4 _ e\2,,2 



1 + A^r „ 8^2^2 



-Sky — 2J||7k 



Bkz = 2J||7k 



(1 + A2)' (1+^2)2 



(1-A^^)+A2(1-m2)- 
(l + A2)(l+/i2)2 ' 

^1 + A2) + (1 - A2) ^4 + 2^1^ 

(1 + A2) (1+^2)2 



9l = b-' /i' 1 



1+^2 ^1^A2 



2, ,3 



8J,i 



A2/i 



i 



= -4^||7k 



(1 + A2f (1 + /X2)- 



(l-7k) 



(l + A2f/2(l + ^2) 



1 + A' 



1+Ai2 



In the zero field case {B = 0, fi = Q) we have Ck = -Dk = 
0, and therefore flkx = i^kx, ^ky = ujky- In the FPF 
phase {B > Bc2, \ — oo, fi = 1) we find i?k = Ck = 
Dk = 0, and the Hamiltonian (|l^ ) is already diagonal in 
the ika operators. 
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